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Abstract 



Search for High Energy Neutrino Emission from Gamma-Ray Bursts with the 
Antarctic Muon and Neutrino Detector Array (AMANDA) 

by 

Ryan Carlton Bay 
Doctor of Philosophy in Physics 

University of California at Berkeley 
Professor P. Buford Price, Chair 



The photo-meson production of pions by shock-accelerated protons could generate a burst 
of ^ 10^^ eV neutrinos from gamma-ray bursts (GRBs) observable in the Antarctic Muon 
and Neutrino Detector Array (AMANDA) or its larger successors. Measurement of this 
flux can test the hypothesis that GRBs are the sources of the highest-energy cosmic rays, 
and GRB neutrinos could permit high-precision experiments in neutrino limiting speed, 
neutrino oscillations, and the weak equivalence principle. Neutrino emission can be expected 
primarily during the prompt gamma-ray flash and satellite coincidence provides a well- 
defined window in position and time that can be searched for an excess of upgoing muon 
events in AMANDA from bursts in the Northern Hemisphere. Using an event quality 
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analysis to further reduce background in a sample of 78 GRBs from the 1997 AMANDA- 
BIO data set, I find a fluence limit of < 3.8 x 10-^min{l, £;^/£;break} TeV cm"^ 
per average burst, which is orders of magnitude more stringent than in similar previous 
searches. 



Professor P. Buford Price 
Dissertation Committee Chair 
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Chapter 1 

Gamma-ray bursts 

...and heaven will bless those stars 
that shine not the longest, 
but those that shine the brightest. 

Anonymous 

The spacecraft Vela 5B was placed roughly 180 degrees apart from its twin, Vela 
5A, at a geocentric distance of 118,000 km on May 23, 1969. The US satellites were part 
of a classified series equipped to aid in enforcement of the Partial Test Ban Treaty of 1963. 
Any nuclear test detonations in the upper atmosphere or on the dark side of the moon 
would have generated signals in the scintillation X-ray, gamma-ray, and neutron detectors 
on-board. Although no clandestine explosions were ever confirmed, during the ten-year 
period April 1969 to April 1979, the Vela group observed 73 instances of a very puzzling 
new phenomenon. Flashes of localized gamma radiation from seemingly random directions 
would suddenly outshine all other sources of gamma rays combined for a few seconds and 
then disappear. The discovery of gamma-ray bursts (GRBs) was announced in 1973 Q. 

The new phenomena generated tremendous excitement over the next two decades 
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and theories of GRB origin flourished. Some placed them at cosmological distances, others 
involved Galactic neutron stars, and still others reduced them to mere twinkles from the 
Oort cloud of our solar system. The picture did not improve significantly until the launch of 
the Burst and Transient Source Experiment (BATSE) on-board the Compton Gamma-Ray 
Observatory in the spring of 1991. BATSE confirmed the earlier suspicions that GRBs are 
very isotropic on the sky and also revealed a deficiency of faint bursts. These facts eliminated 
the possibility that GRBs originate on neutron stars in the Galactic diskQ, and made the 
case for a cosmological distance scale compelling [§, ^J. The science of GRBs was recently 
revolutionized and the cosmological hypothesis confirmed with the launch of the Italian- 
Dutch BeppoSAX satellite, which allowed counterpart observations to localize several GRBs 
to host galaxies and measure emission and absorption line redshifts of z ~ 1 jQ, ^ |5|. 

1.1 GRB phenomenology 

If GRBs emit their radiation isotropically, then their cosmological distance scale 
implies that they are extremely powerful and rare events, releasing in a matter of sec- 
onds 10^1 to lO^-^ ergs in gamma rays alone. Recent observations of afterglow light curve 
breaking 0, ^ suggest beaming solid angles of order ^7 0.1, which would reduce these 
energies by factors of (|^)~^ ~ several hundred. Assuming a constant formation rate, the 
observation of ~ 1 per day at Earth implies a rate of only 10~'^(|^)^^ yr^^ galaxy"^. 

One reason GRBs provoke such wild speculation is the enormous variety they 

exhibit. In so many words, they are short flashes of low-energy gamma rays with a non- 
^The BATSE result strongly constrained but did not rule out an extended Galactic halo origin. 
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thermal spectrum. Observed durations, however, span at least five orders of magnitude. 
Light curves are complicated and vary drastically from burst to burst. The spectrum is 
very hard, with emission out to GeV Q and even TeV |ll| energies. There is also now 
an indication of two distinct classes of GRB, with appreciably different distance scales. 

1.1.1 Spectra and light curves 

The GRB photon spectrum can generally be well fit 
joined smoothly at a break energy of ~ 200 keV: 

{{hu)'^e^~^^ for hv < H 

[(a-/5)^o]^°"^nH^e(^""^ for hv>H, 
where H = (a — f3)EQ. GRBs typically emit in the few hundred keV energy range, while 
prompt emission in the X-ray and lower energies is very weak. The spectrum is exceptionally 
hard, and most bursts have a high-energy tail with E'^N{E) nearly constant. On the low- 
energy side the spectrum often follows Fj^ocz/",— ^ < a < |, which is consistent with the 
tail of synchrotron emission from relativistic electrons. The prompt gamma radiation is 
followed by x-ray, then optical, then radio emission over a time scale of days. 

The shortest recorded GRB lasted a mere 5 milliseconds. The longest, GRB940217, 
showed GeV activity 90 minutes after the initial outburst. GRBs 961027a, b, c, and d, the 
so-called "gang of four", occurred over a period of 2 days with overlapping error boxes. If 
these are in fact the same object, then GRB durations span an astonishing seven orders of 
magnitude around an average of ~ 10 seconds. 

The distribution of burst durations is bimodal with peaks at 0.3 and 25 seconds. 
Bursts can be clearly divided into those with duration greater than and those less than 2 



12 1 with a double power law 



seconds, with a now definite tendency for the shorter subset to have harder spectra. There 
is a clear correlation between burst duration and hardness^ for the entire set of bursts, while 
the two quantities are not at all correlated within the two subsets |l3| . 

Most GRBs show time structure 5T <ti T, their total duration. The shortest rise 
times recorded are roughly equal to the shortest structures within time histories, ST^[^ ~ 0.2 



msec |14]. This is one of our most important observational facts, since it constrains the size 



of the emitting region to the scale of a compact object: 

-Remission < cdT^ain ~ 60 km 

1.1.2 Progenitors 

All of the most fashionable theories of GRB origin eventually culminate in the 
formation of a black hole surrounded by a torus of debris, in order to take advantage of the 
potentially short time scales and large extractable energy reserve, ~ 10^^ ergs, which such a 
system provides. Afterglow observations of some bursts suggest dense gaseous environments 
characteristic of star-forming regions, giving rise to "hypernova" models |]l5| involving core 
collapse of a massive star, as explanations for the longer and softer bursts. Mergers of 
binary neutron stars or neutron star-black hole systems popularly used to explain shorter 
GRBs are expected to involve older stellar populations in more rarefied environments, the 
violent formation of the binaries having propelled them to distances of several kiloparsecs 
from their birthplaces. Compact binary mergers would be expected at redshifts 50% — 80% 
those of hypernovae. 



^GRB hardness is most often defined as the fluence [ergs cm ^] in the fOO — 300 keV channel divided by 
the fluence in the 50 — fOO keV channeL 



Regardless of how the black hole-torus system forms or the mechanism through 
which its gravitational or spin energy is tapped, any respectable model must be able to 
produce sensational energy release. When combined with observed luminosities, the size 
constraint above implies a region with an enormous photon energy density of about 10^^ 
erg cm~'^. Such a plasma has very large optical depth to pair production, yet observation 



requires that it be optically thin. It was first pointed out in the 1980s by Goodman |1£] and 
Paczyhski |]l^] that this problem can be solved if the emitting region is moving relativistically 
toward the observer with a sufficiently large Lorentz factor T. Then the photons at the 
source are of energy E'observed/rj and the size constraint is relieved by a factor F^: 

-^^emission ^ T c6Tmin- 

Infeasibly high energies are required for the entire source region to move relativistically, 
solely to overcome this "compactness problem." The most economical scenario, both con- 
ceptually and energetically, is one in which the relativistic motion is intimately related to 
the observed burst. This is achieved naturally by the relativistic fireball shock model. 

1.2 Great balls of fire 



In the model first developed |2| by Meszaros and Rees [|^, Narayan et al. ||19| 



and Paczyhski and Xu |20|, GRBs are the result of deceleration of ultra-relativistic matter. 
The model is essentially three-stage: 

i) An inner engine produces a highly variable relativistic outflow of (optically thick) 
plasma. 

ii) This energy is transported via bulk motion of the plasma out to ^ 10^^ cm before 
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it becomes optically thin. 

iii) The energy is dissipated, emitting the prompt gamma rays and later radiation at 
longer wavelengths. 

The plasma acceleration must saturate at a Lorentz factor Fmax ~ 'm~~ — ^ of least 100 
for the gamma rays to escape pair-production losses, which constrains the mass of the flow: 

Moutflow < IO-^Mq. (1.1) 

All of the energy of the fireball will be converted into bulk kinetic energy of any con- 
tamination of baryons entrained in turbulent magnetic fields. Energy conversion is then 
accomplished through shocks and the observed photons are the synchrotron or synchrotron 
self-Compton radiation of relativistic electrons which have undergone first-order Fermi- 
acceleration. 

The prompt gamma radiation is most likely due to internal shocks in which faster 
shells of the outfiow overtake slower ones. The progenitor cannot produce a single explosion 
but rather must generate a wind, while the burst duration and variability are determined 
by those of this inner engine. 

External shocks also occur when the blast wave encounters the surrounding medium. 
It is possible to explain the rapid variability seen in the gammas if the density of the sur- 
rounding medium is highly irregular, but this is prohibitively inefficient. There is also no 
evidence that the width of features in the gamma-ray light curves increases with arrival 
time as the fireball decelerates. External shocks are instead probably responsible for the 
afterglow at longer wavelengths over time scales of days and weeks. 
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1.3 Future missions 

Our understanding of the nature of GRBs is certain to improve over the next 
several years as a new generation of satelhtes permits not only better measurements of the 
initial outburst, but also better measurements of the afterglow by an increasing number 
of ground-based observers. Currently GRB afterglow observations typically begin >> 10 
seconds after the gamma-ray flash. Fast multi- wavelength follow-up observations are needed 
to constrain parameters of the fireball and of the surrounding medium. The two most 
notable future missions, HETE-II (High Energy Transient Explorer) and Swift, promise to 
provide GRB analysis in almost real time. HETE-II will be launched in July 2000 and will 
quickly transmit arcminute error boxes for about 50 bursts per year. Swift is expected to 
launch in 2003 and will be equipped with a wide-field X/gamma-ray camera, a focusing soft 
X-ray telescope and an optical telescope. It will quickly slew to the positions of about 300 
bursts per year for detailed X-ray and optical work [|l4| . 
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Chapter 2 



Neutrino astronomy 



Is to the sky 
Looking to the sky 
and down 
Searching for a ground 
With my good eye closed. 

Soundgarden, Searching with my good eye closed 



2.1 The neutrino as an experimental device 

Information about distant objects is delivered to us by light from radio frequencies 
up to the highest measured gamma rays provided that the photons have a clear path. 
We cannot view directly radiation originating within astrophysical objects or from regions 
obscured by dust. At high energies, gamma rays suffer pair-production losses on other 
photons whenever 

^-Egamma^photon ^ (^JTlelectronC ) ■ 

This means that PeV gamma-rays do not make it more than a few kiloparsecs before being 
absorbed by the cosmic microwave background radiation (CMBR) , and TeV gammas could 
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be absorbed within 100 megaparsecs by diffuse infrared. 

Cosmic rays provide vital information as well, but are not very useful for a tele- 
scope. A proton with energy less than about 10^^ eV does not point back to its source 
because its direction has been randomized by intervening magnetic fields (the field of 
our Galaxy is enough). Above this energy it will interact with the CMBR through the 
A-resonance, since 

4ii/proton £photon ^ ('^A '^proton • 

Neutrinos could serve as a useful diagnostic of interesting places in the universe. 
The high-energy neutrinos from the decay of mesons will trace the spectrum of those cosmic 
rays interacting either with interstellar gas or possibly even at their acceleration sites. High- 
energy neutrino production is also expected in the decay or annihilation of ultra-heavy 
particles. Such particles might be concentrated at the core of the Earth, Sun, or Galactic 
center, where relic weakly interacting massive particles (WIMPs) produced during the Big 
Bang have fallen into a gravity well and been trapped via elastic scattering. 

The neutrino can also escape its birthplace and traverse vast distances unimpeded. 
With its tiny mass and inability to take part in any but the weak interactions, the interac- 



tion length of a 1 TeV neutrino, for example, is about 2.5 million kilometers of water |21|. 
Unfortunately these same properties that make the neutrino an immutable courier of in- 
formation make it diabolically elusive to detection. Nonetheless, the neutrino telescope 
can look to several experimental precedents for encouragement. Solar neutrinos have been 



observed both radiochemically |22] and through the observation of recoiling electrons |23|. 



Low-energy neutrino astronomy was christened with the detection of neutrinos from super- 
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nova SN1987A |2^, |25|. More recently, particle physics using atmospheric neutrinos has 
become a serious science, with convincing evidence for neutrino oscillations [ p^ . 

2.2 Detection method 

With cross-sections around lO"(^^^^^y"^°™^^^™§) cm^, any experiment bent on 
neutrino detection must patiently await and somehow discern the particle's rare interactions. 
In one conspicuous reaction of great interest for the purposes of a neutrino telescope, the 
Ufj^ produces a muon through the charged-current interaction. This daughter muon carries 
away one-third to one-half of the neutrino energy, and at high energies preserves the neutrino 
direction with a mean deviation of 

^(6*2^) [radians] 

which is about one degree at 3 TeV [^^. This muon can be traced by the ~ 250 Cerenkov 
photons it emits per centimeter as it traverses a transparent medium embedded with a 
lattice of photomultiplier tubes (PMTs) [^^. Unlike other much smaller water Cerenkov 
detectors which have typically a factor of 10^ higher density of PMTs, the neutrino telescope 
is designed primarily for the TeV range, taking advantage of the long lever arm furnished 
by a high-energy muon track. A muon in ice or water loses only about 2 MeV cm^^ to 
ionization up to 5 TeV, above which the energy loss rises approximately linearly with energy 
due to pair production, bremsstrahlung, and nuclear interactions. A muon which does not 
suffer catastrophic energy loss can reach the instrumented volume of a detector several 
kilometers away from a high-energy neutrino interaction within a large effective volume. 
The electromagnetic showers induced along the tracks of the higher-energy muons add to 
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their light generation, increasing detector effective area with energy as weh. Calorimetry is 
possible albeit imprecise, and depends strongly on PMT spacing and the optical properties 
of the medium. 




Figure 2.1: {Left) A relativistic muon can be traced by its Cerenkov cone as it traverses a 
lattice of photomultiplier tubes. (Right) An electromagnetic or hadronic shower event. 



A lattice of PMTs provides other modes of neutrino detection. Reconstruction of 
the energy and perhaps even direction of contained cascade events arising from electron neu- 
trinos, stopping muons and the "double bang" signature of a tau neutrino is possible. The 
low-energy neutrino burst from a nearby supernova can be noticed by carefully monitoring 
only the counting rate of the PMTs |28]. 

At low energy thresholds a neutrino telescope must cope with the comparatively 
enormous flux of muons produced in the atmosphere overhead. The simplest strategy is to 
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accept only upward traveling muons: an upgoing muon is certain to have been a neutrino 
which penetrated the filter of the Earth. The detector must then be able to distinguish 
an upward- from downward-traveling event very reliably. An above-ground detector would 
need better than 10^^ discrimination at sea level, while a detector shielded under 1000 
meters of water still needs about 10^. 



Figure 2.2: The neutrino telescope must determine muon direction reliably in order to 
distinguish extraterrestrial neutrinos from atmospheric backgrounds. 

Neutrinos are also produced in cascades from particle collisions in the atmosphere 
on the other side of the Earth. In general the flux of neutrinos from pion decay is given 



where At^i, = Znt^{1 — r.,^)'^ /{'^ + 1), r^^ = {m^/m.,^)'^, Bj^^ is a constant that depends 



extraterrestrial 
neutrinos 




atmospheric 
|/ neutrinos 



by© 
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on the spectral index 7 and nucleon and pion attenuation lengths, and the (• • •) term 
represents the contribution from other mesons. The energy e^r is a characteristic quantity 
that reflects the competition between decay and interaction in the medium, and for cascade 
development in the Earth's atmosphere ~ 115 GeV. Although there is no direct way 
to distinguish these from the extraterrestrial variety, atmospheric neutrinos have a smooth 
angular distribution and above several TeV their spectrum falls precipitously as E~'^'^, at 
least out to an energy above which charm production (with a flatter spectrum) dominates. 
In investigations of localized transient phenomena, background from atmospheric neutrinos 
is virtually negligible. 




log, „ Energy 



Figure 2.3: Schematic presentation |27] of neutrino fluxes from (1) atmospheric neutrinos, 
(2) cosmic-ray interactions with interstellar gas and (3) cosmic-ray interactions at their 
source. 



The + A ^ fi + X cross-section also grows linearly with energy up to about 
10^^ eV and then levels off to increase as E^''^ when the mass of the W propagator begins 
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to impede its coupling. Neutrino absorption in the Earth increases until eventually looking 
downward through the planet becomes a disadvantage. Above about 40 TeV, the neutrino 
interaction length becomes smaller than an Earth diameter, about 10^ m of water equivalent. 
The stiffer the neutrino spectrum, the more events are confined to solid angle closest to the 
horizon. 
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Chapter 3 



AMANDA 



The Antarctic Muon and Neutrino Detector Array is an American and European 



effort to construct a neutrino telescope at Soutli Pole Station using the Antarctic ice cap |29|. 
This ice is a stable platform for the delicate instruments of a particle experiment, and 
also an extremely transparent optical medium with manageable scattering for faint light 
detection with sparsely distributed PMTs. The unique location of AMANDA permits (and 
in some cases requires) study in other scientific disciplines including glaciology, seismology, 
paleoclimatology, and even microbiology. 

Construction of the 10-string phase of the array, AMANDA-BIO, was completed in 
the austral summer 1996-97 and its geometry and calibration established in 1998. AMANDA- 
BIO has a threshold of about 50 GeV and achieves an angular resolution of 2.5 degrees, while 
presenting an effective area on the order of 10, 000 m^ to the highest-energy neutrinos. The 
data set from the 1997 season of this first-stage prototype will be the primary subject of 
this analysis and is relatively well understood. Significant progress into calibration of the 
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1998 and 1999 data has been made. 



3.1 Strings of beads 

The BIO array consists of 302 PMTs arranged on 10 strings at depths between 
1.5 and 2 km in approximately a right cyhnder of 60 meter radius, a geometry biased to 
vertically-travehng muons. The acceptance of the array was made more symmetric with 
the addition of three strings in 1997-98 and six more strings in 1999-00 for a total of 19 
strings and a 90 meter radius. This AMANDA-II will comprise the core of the IceCube 



kilometer-scale array still in the proposal stage at the time of this writing [30|. 

The maximum construction rate of the array achieved thus far has been 6 strings 
per season. A large hot-water drill is used to make a 60 cm diameter wet hole down to a 
depth of 2 km, and the AMANDA string is then deployed before the hole refreezes. Complete 
refreeze typically takes several days due to residual heat. An optical module (OM) consists 
of an 8-inch PMT and base within a high-pressure glass sphere, with its own electrical 
connection to the surface for power and signal transmission, and a fiber optic connection 
for calibration purposes. Signal processing is carried out entirely on the surface in the 10- 
string prototype. Although this configuration is quite robust and the low background rates 
of PMTs in ice allow for a very simple triggering scheme, the collaboration is exploring many 
more sophisticated technologies for use in the IceCube array. Fiber optic signal transmission 
and full waveform digitization can improve photoelectron timing and dynamic range. Local 
coincidence triggering between neighboring OMs might become necessary in place of the 



simple multiplicity trigger as the number of OMs continues to grow |31|. 
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ice surface level 




Figure 3.1: The first four strings were deployed too shallow in ice contaminated by air bub- 
bles which have not yet converted into clathrates (air-hydrate crystals), and these comprise 
AMANDA-A. AMANDA-B was deployed much deeper in bubble-free ice. An individual 
optical module is exploded to illustrate its mechanical design. 
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3.2 An all-natural optical medium 

The southern polar ice cap is the product of hundreds of thousands of years of 
steady accumulation of less than a foot per year of very finely powdered snow. Ground- 
penetrating radar isochron maps taken from airplanes suggest that, on the kilometer scale, 
South Pole ice can primarily be described as a function of depth. Layers of higher con- 
centrations of various kinds of dust correspond to glacial maxima (ice ages), when the 
atmosphere was more arid and the larger temperature gradient between the poles and the 
equator created higher winds. 

AMANDA has made possible the study of the optical properties of South Pole ice 



at depths from 800 to 2300 meters [32, 33, p3, B^, 3(:, 37]. Isotropized and filtered light from 



pulsed LEDs and lasers, DC LEDs, halogen lamps, and xenon arc lamps deployed alongside 
or inside the optical modules, and laser pulses transmitted via fiber optic from the surface 
serve as beacons to OMs as far as 200 meters away. 

A diffusive analytic solution can be applied effectively only when the source and 
receiver are separated by a distance d much larger than one effective scattering length, Ac = 
Ascat/(1 — (cos 6)). If the scattering length is much larger than the spacing (the ideal case for 
a telescope), delay time is a delta function. Monte Carlo simulation becomes necessary to 
describe bulk properties when d ~ Ac, as is the case at depths > 1400 meters in AMANDA. 
Pulsed sources can be used to extract separately both the effective scattering length and 
absorption length, Aa, from arrival-time distributions by fit to simulation. The propagation 
length Aprop = \/ AaAc/3 can be determined over large distances from a measurement of 
fluence versus distance falling as N{d) oc 4e ^p"^"? . 
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Emitter/receiver combinations in the same horizontal plane at as many depths 
as possible are desirable to avoid integrating over vertically varying features such as dust 
bands. Optical properties can be predicted all the way down to bedrock by matching 
variations measured with optical methods to dust concentration profiles of ice cores taken 
at the Russian Vostok site ~ 1300 km away and the Japanese Dome Fuji site ~ 1500 
km away and extrapolating (Figure \i.2\ ) ||3^ . More is now understood about the optical 
properties of the Antarctic ice cap than is known about those of the deep ocean Is^. 



Figure |3.3| is a plot of South Pole ice absorption. Shown are data fitted with a 
"three-component" model of Bergstrom and Price ||3^ . Absorption increases dramatically at 
wavelengths shorter than the visible due to electronic excitation modes and at wavelengths 
longer than the visible due to molecular modes. Inside this window the transparency of 
South Pole ice is limited only by the small concentrations of dust, with an absorption 
length on the order of 100 meters. This is many times more transparent than laboratory ice, 
and combined with the moderate concentration of scattering impurities which act as wave- 
guides. South Pole ice provides a natural medium for efficient calorimetry of electromagnetic 
processes. 
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Figure 3.2: Matching the scattering coefficient measured opticaUy at South Pole to ice 
core data at Vostok station determines the relation between age and depth. LGM is Last 
Glacial Maximum, while D signifies a lesser cold period not generally considered an ice age. 
A, B, and C are still smaller anomalies used as reference points. 
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— Bergstrom-Price 3-component model 




wavelength [nm] 



Figure 3.3: The three-component model of South Pole ice absorption. In the visible window 
between 4 x 10^^ Hz and 8 x 10^^ Hz the transparency is limited only by a small level of 
dust contamination. 



The simulation [40| in Figure 3.4 demonstrates how the response of a Cerenkov 
detector can be affected when optical properties are not uniform. Layers of enhanced concen- 
trations of various impurities can either decrease or increase the light-collection efficiency of 
the array. The effective scattering length is found to vary by a factor of two over AMANDA- 
BIO, yet all currently distributed collaboration software (including all software used in this 
analysis) assumes uniform optical properties for the sake of simplicity and speed. This is 
a primary cause of current disagreement between data and Monte Carlo simulations, and 
introduces considerable uncertainty into any estimation of sensitivity. 
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Figure 3.4: ( Top) Monte Carlo simulation in which optical properties are fitted to an 
arrival time distribution between an emitter and receiver separated by 50 meters at varying 
distance from a dust layer with three times more scattering and absorption. Near the 
layer, optical properties are actually perceived to improve, due to the layer acting as a 
backboard. {Middle) Monte Carlo simulation of the response of an array of optical modules 
to downgoing muons when an impurity layer is introduced with three times more absorption 
and scattering and (bottom) with increased scattering only. The dashed lines show response 
with uniform ice. In the scattering-only case the layer acts as a photon guide and enhances 
the hit rate on the upstream side. 
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3.3 Event reconstruction 



The arrival time distribution of a population of photons emitted from a pulsed 
point source is a combination of delay due to scattering impurities and the usual Gaussian 
timing uncertainties of several nanoseconds due to electronics (Figure 3.5). Because a neu- 
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Figure 3.5: The arrival time distribution of a pulsed blue LED at a distance of 100 meters. 



trino telescope is sparsely instrumented and AMANDA-BIO is still small in comparison to 
the range of the muons it is designed to detect, a likelihood analysis is necessary instead 
of a simple minimization to make the most of degraded information delivered by de- 
layed photons. The quantity Pi{t) gives the probability of OM i being hit at time t, and 
reconstruction seeks to maximize |41] 



l0g(/:) = logdlallhitsl'i) = Eallhitslog(pi). 

About 5% of downgoing muons are mistaken for upgoing by the reconstruction, 
and quality criteria are necessary to identify these events automatically. An effective cut 
parameter is Ni^ij.(,(.t, the number of OM responses presumed to be minimally scattered pho- 
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tons from the Cerenkov cone of a given postulated track, to a specified degree of directness. 
An event with many direct hits is one with a comparatively high content of uncorrupted 
information. The number of direct hits is classified by A'dij-A, -^dirB) and A'^^irC to distin- 
guish OM responses which occurred within [—15,25], [—15,50], and [—15,75] nanoseconds 
of the expected time. Combining these with various topological cuts improves up-down 



discrimination (Figure p.6D and angular resolution to acceptable levels. Figure |3.7| shows 
the angular distribution of the remaining neutrino candidates in the 1997 data passing a 



full set of quality cuts |4^ . 




0.6 0.8 1 
cosine(zenith) 



Figure 3.6: Reconstructed zenith angle at three cut levels on the same sample of events. 
Requiring a minimum number of direct hits or, e.g., a minimum span of hits along a track 
(■Zhits) removes events misreconstructed as upgoing. 
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Figure 3.7: Skyplot of 1997 BIO neutrino candidates from one analysis in Galactic coordi- 
nates. The detector has a significant vertical bias which is further exaggerated by quality 
criteria. 
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Chapter 4 



Neutrinos from hell 



/ am fire and air; my other elements I leave for the baser forms of life. 

William Shakespeare, Anthony and Cleopatra 

Gamma-ray bursts can be regarded as the relativistic analogue of supernova rem- 
nants (Table |4.lD . A highly efficient, relativistic outflow produces gamma radiation in 
internal shocks in the case of GRBs, whereas the non-relativistic outflow from supernovae 
produces optical emission in external shocks with notoriously low efficiency. The analogy 
can be taken a step further: just as supernovae are invoked to explain cosmic-ray energies 
up to 10^^ eV and even beyond, GRBs represent a prime candidate location for cosmic-ray 

Table 4.1: GRBs are the relativistic analogue of supernova remnants. 



SNR 


GRB 


high 


Optical depth 

low 


Outflow velocity 
non-relativistic (F < 1.02) relativistic (F > 100) 


low 


Efficiency 

high 




Emission 


optical 


gamma rays 




Shocks 


external 


internal 
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acceleration to energies beyond the "ankle," the flattening in their spectrum near 10^^ eV. 
A compelling fact is that the inferred energy injection in gammas by GRBs, about 10^^ erg 
Mpc~'^ yr~^) is strikingly similar to that needed to produce the ultra-high energy cosmic 
rays (UHECRs) ||. 



4.1 Hadronic GRBs 



The mass limit in Equation is often referred to as the "Baryon Contamination" 
problem because it is difficult to develop a scenario in which only tiny quantities of material 
are entrained in the outflow of the fireball. A small amount of baryons is actually useful 
in explaining how the fireball evolves from radiation-dominated to matter-dominated before 
becoming optically thin (Section |1.2[ ). 

In portions of the flow where electrons are accelerated, protons could undergo 
shock acceleration as well, possibly to ultra-high energies. Because the particles spend so 
little time on the shock front, it is difficult to efficiently accelerate protons on the highly 
relativistic external shocks. Internal shocks, on the other hand, are expected to be only 
mildly relativistic, since Lorentz factors of different shells should differ by factors of order 
unity and not more than an order of magnitude. If a sufficient fraction of the wind energy is 
carried by magnetic fields which can isotropize proton velocities on either side of the front, 
and if the proton acceleration time is shorter than both its synchrotron loss time and the 



wind expansion time, GRB fireballs could accelerate protons to > 10^'' eV p4|, I45[. There is 



now evidence for TeV [IC] and sub-TeV [11 1 gamma emission from GRBs which cannot be 
explained by a simple extrapolation of the BATSE energy spectra, but could be the result 
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of ttq production or synchrotron emission by energetic protons |T^]. 

4.2 The z/ flux 

Protons accelerated in GRBs will photo-produce pions in interactions with fireball 



photons and the decay of these pions could generate a burst of high-energy neutrinos 47 1: 



p + 7 — > A+ 
\ 

vr"*" + n 
\ 

\ 

The predicted energy distribution of particles accelerated on strong shocks is a power law 
with Np{Ep)dEp oc Ep ^'^^''^dEp. These will interact with GRB photons to produce neutri- 
nos which trace the broken power law of the seed photons: 

N„{E,y)dE„ oc — dE^, with /3 = 

-'-'break 



and 



2 for El, > E'break 



^break ~ 7 X lO^^Tioo eV, (4.1) 



where Fsoo assumes a fireball Lorentz factor of 300. 

Most calculations take advantage of the A-approximation, which assumes the 
photo-production cross-section is dominated by the A-resonance near threshold. At high 
energies, non-resonant interactions of protons with the flat part of the photon spectrum 
could become significant [|^. Any neutrino emission is probably suppressed above about 



29 

10^^ eV anyway, however, when the hfetimes of the intermediary muons and pions within 
the cascade exceed their synchrotron and adiabatic loss times. 

The neutrino flux is determined by the factor /ttP, where p is the injection rate of 
high-energy fireball protons and /jr is the fraction of their energy given to pion production. 
If GRBs are responsible for the UHECRs this fixes p ~ 10'*'^ erg Mpc^'^ yr^^. Taking 
/tt = 0.15 the energy of GRB neutrinos can be normalized, 

by determining the constant A, found by Halzen and Hooper |49| to be 1.2 x 10^^^ TeV 
cm~^ s~^ sr~^. Assuming 1000 bursts per year over 47r steridians, the fluence of neutrinos 
from a single GRB, in the compact form of Waxman and Bahcall, is then 

El^ = 4.8 X 10- Vin{l, ^./^break} TeV cm'^. (4.2) 

This is consistent with the neutrino fluence estimated directly from the gamma rays [^ ]: 

ii;2^ = 8.5xlO-^(:^)TeVcm-2, 

where /e is the ratio of the energies of the electromagnetic and hadronic components of the 
fireball. 

Event rates in a neutrino detector can then be estimated by a convolution of 
the flux and detection probability of muon- flavored neutrinos, Pu^jj.-, determined by the 
interaction cross-section and muon range: 

This predicts that a detector with 1 km^ effective area should see on the order of 10 events 
per year, which is discouragingly small. 
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Table 4.2: A values for various models. In the Halzen/Hooper and Waxman/Bahcall 
predictions burst parameters are fixed at typical values. The model of Halzen, Alvarez and 
Hooper includes fluctuations in energy and distance. 



Model 


A 


-^brcak 


Halzen and Hooper 
r = 300 


1.2 X 10"^^ 


700 TeV 


Waxman and Bahcall 
r = 300 


4.5 X 10"^^ 


500 TeV 


Halzen, Alvarez and Hooper 
r = 300 


6.2 X 10-^^ 


700 TeV 


Halzen, Alvarez and Hooper 
r = 100 


9.1 X 10"^^ 


77 TeV 



It is heartening to note that this flux is the prediction for the typical GRB. It may 
be far more likely to detect neutrino emission from a single burst with highly favorable 
properties than many standard-candle bursts. An exceptionally close and powerful burst 
could of course provide better-than-average observational possibilities. Fireball parameters 
are not well constrained and the photon target density and hence the value of /jr and 
the break energy all depend sensitively on the Lorentz factor F of the outflow. Burst-to- 
burst fluctuations could increase the neutrino intensity from a particular burst, perhaps 



dramatically higher than that in Equation 4^ |4£, 51 1. Some eff'ective A values, which take 



into account fluctuations in energy and distance, are shown in Table 4.2. Integrating the 
flux of GRB neutrinos from ii^min, the 1.232 GeV threshold of the A-resonance, the fluence 
is then proportional to 

^2 ^ A -E'break _^ 

dEu -E'break -^min 

More detailed observation will constrain the degree to which parameters of the GRB flreball 



fluctuate and how those fluctuations could be manifest in any neutrino emission. 
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4.2.1 Scientific goals 

A measurement of the GRB neutrino flux can test the hypothesis of GRBs as the 
origin of the highest energy cosmic rays. Like SN1987A, GRBs could permit even stronger 
tests of the neutrino limiting speed and the weak equivalence principle at levels of 10"^^ and 
10~^, respectively, due to their short duration. Furthermore, since no v-r- are expected to be 
produced in the source, an appearance of would be clear evidence of neutrino oscillation 
at a value of Am^ as low as 10^^^ eV^. 

4.2.2 Monte Carlo simulation 



The stiff GRB neutrino flux should suffer significant Earth shadowing |21|. Monte 



Carlo simulation |52] gives the expected event rate at three cut levels (Section |5.2| ) in 
AMANDA-BIO versus zenith angle for 500 Halzen and Hooper GRBs distributed over 27r 



steridians in Figure |4.1| . The sensitivity of the array drops by nearly an order of magnitude 



near nadir, where effective area and angular resolution are otherwise optimum. 
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Figure 4.1: Monte Carlo calculation gives expected event rates from 500 GRBs spread over 
10 zenith bins. The dependence is dominated by neutrino absorption in the Earth. 



33 



Chapter 5 

Data processing 

5.1 Data collection 

5.1.1 GRB data 

Localized burst data are available from three sources, and I compiled triggers in 
the following priority: 

i) BeppoSAX triggers, 

ii) bursts detected by three or more Inter-Planetary Network (IPN) satellites, 

and 

iii) BATSE triggers. 

The BeppoSAX Wide Field Camera covers about 5% of the sky and together with the 
Narrow Field Instrument provides arcminute locations to a few bursts per year. Bursts 
detected by satellites in the IPN must be observed by at least three widely separated 
spacecraft for their position to be well localized. While there are many spacecraft near 
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Earth, in interplanetary space there were only two in 1997, Ulysses and the Near Earth 
Asteroid Rendezvous (NEAR). BATSE remains the primary source for GRB triggers as of 
this writing BATSE triggers are based mostly on counts in the 50 to 300 keV range and 
it is very possible that a large population of hard GRBs escapes detection. Searches of 
BATSE data have produced many hundreds of non-triggered burst candidate events |]53|| , 
most localized to no better than 10°. Some events have been confirmed by the Ulysses 
experiment, while some others are probably spurious [|l^ ]. 

BATSE consists of eight modules of Nal(Tl) scintillation detectors, one module on 
each corner of the Compton-GRO satellite. BATSE can estimate source direction by the 
relative intensities in the eight modules and the BATSE team reports the statistical error 
radius of a 68% confidence ellipse. For a weak burst, the positional error can be as large as 
20°, but for strong bursts it is less than one degree. There is in addition a systematic error 
of about 1.6°(73% probability) [^]. The BATSE team uses Tqo{T5o), the time required 
to accumulate 5%(25%) to 95%(75%) of the total counts, to characterize burst duration. 
Those bursts which suffer data gaps during the event (about 15% of the triggers in 1997) 
are not assigned TgoiT^o) times, although some visual estimates of duration are available in 



the BATSE Comments Table g. 



5.1.2 AMANDA data 

I extracted all two-hour intervals of AMANDA-BIO muon data bracketing 1997 
GRBs from the collaboration repository after the Level 1 filter. These events have been 



^BATSE and the other instruments on-board the Compton Gamma-Ray Observatory were turned off at 
UT midnight on my 28*'' birthday, May 26, 2000, due to fears of loss of spacecraft gyro-control. Compton- 
GRO was removed from orbit and burned up in the Earth's atmosphere on June 4, 2000. 
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calibrated and reconstructed with both a simple line fit and the full reconstruction, while 
ignoring malfunctioning optical modules and isolated hits presumed to be noise. The line 
fit result is required to have 9 > 50° for an event to be part of the 22% passing this stage 



of the filter H]. 



I first analyzed these data to remove periods of uncharacteristic detector perfor- 
mance from the search, isolating hardware and software glitches by consulting previous data 
diagnostics and using a specifically designed algorithm. Each two-hour measurement 
was divided into ten-second intervals and evaluated on the basis of total array rate and 
occurrence of long gaps between events. Time bins in which the total array rate deviated 
by more than 3.5 to 4 standard deviations from a running average or which contained more 
than two gaps longer than 0.5 seconds were removed. Bins immediately next to a failing 
bin were required to pass even stronger criteria. The cuts were tuned by visual inspec- 



tion to minimize biasing. Figure 5.1 shows an example application of the algorithm and a 



calculation of live time. During a few bursts the muon data were completely unusable. 



Figure 5.2 shows how the background of noise and misreconstructed downgoing 
events is removed as more direct hits are required in the events. The cut A'^^irB ^ 2 results 
in roughly the same background suppression as A'dirC ^ 3, and A'^dirB > 4 is comparable 
to A^dirC > 7, yet the A^'dirc cut preserves more signal than the A'^dirB cuts in each case 
(Figure O). 



36 




-3000 -2000 -1000 1000 2000 3000 

Level 1 event rate for GRB 6189 in 20 s bins 




0.25 0.5 0.75 1 0.25 0.5 0.75 1 



Time between events (5t) [s] Anomalous 5t [s] 




-3000 



-2000 



loV 



10 

1 



jiud 



-1000 1000 

Excluded periods for GRB 6189 



10 
10 
1 



2000 



3000 



100 200 
Raw rate (20 s bins) 



300 



100 200 
Cleaned rate (20 s bins) 



300 



Figure 5.1: The cleaning algorithm results during a somewhat erratic period of AMANDA 
behavior. Event gaps introduced by the data acquisition system which exceed a fit to the 
natural BIO trigger exponential were integrated to determine live time, which averaged 
about 68% in 1997. Time periods with highly improbable total array rates or excessive 
gaps were excluded from the analysis. 



37 



s. -4 
10 

s 



10 



1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 








































— Ndi 














— Ndi 


yC 








































































































i 






1 1 1 1 


1 1 1 1 


1 1 1 1 


1 1 1 1 


1 1 1 1 


1 1 1 


II 


II 





a 



10 



OS 

o 
a 



6 7 8 9 

Direct Hit Cut Level f >N) 



Figure 5.2: Misreconstructed upgoing background at Level 1, both per degree^ per second 
and in a 10 degree error circle in 10 seconds. 



5.2 Optimization 

Neutrino candidate searches in the current AMANDA typically employ a battery 



of event quality criteria which suppresses signal by 99% [57|. The transient and localized 
nature of GRBs allows quality criteria to be relaxed and the sensitivity greatly increased. 
I optimized cuts using a Monte Carlo simulation of the response of AMANDA-BIO to 
the canonical GRB neutrino spectrum, and AMANDA-BIO data for background. Muon 



propagation was simulated using MUDEDX [jSgl software and the detector response using 
AMASIM p^] . The point spread function of reconstructed GRB events is shown in the top 



part of Figure 5.3. The bottom figure shows how the efficiency increases as the error circle 



is enlarged, with gains in signal diminishing after events accumulate quickly in the first 



several degrees. Figure 5.4 shows how the Q-factor (signal / \/ back ground) and efficiency 
vary together with search bin size. The events were generated from a constant direction 
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near 6 = 45°, cp = 135°. At high energies, angular dependence of detector resolution is not 
expected to be dramatic. 
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Figure 5.3: Top: Monte Carlo simulation of the point spread function of the GRB neutrino 
spectrum in AMANDA-BIO. Bottom: The efficiency as a function of error circle radius at 
various cut levels. 
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Figure 5.4: Top: Q-factor [signal / background) versus search bin size. Bottom: Efficiency 
versus Q-factor. 
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Table 5.1: Selected cuts. 



Ndircct cut 


Bin radius(^) 


Background (10 s) 


Efficiency 


Relative Q (S/y^Bkg) 


None 


31° 


~ 10 


0.61 


1.0 


A'dirC > 5 


22° 


~ 1 


0.45 


2.7 


iVdirC > 7 


12° 


~ 0.1 


0.36 


7.5 



Optimum cuts are those that perform in the uppermost and rightmost parts of 



the graph of efficiency versus Q-factor on the bottom of the Figure 5.4. I chose the three 
cocktails[| in Table 5.1. The tightest cut provides the best discrimination while retaining 
a statistically meaningful number of events and yields the greatest gain when combining 
measurements to improve an upper limit due to the low background, with L cx jp" — as 
opposed to L oc . J- with high background. The Q-factor could be further optimized with 
more sophisticated cuts on Monte Carlo simulations of a better-constrained GRB neutrino 
spectrum, and by making independent zenith and azimuth detector corrections. 



5.3 On-Off source event counting 

I took the satellite resolution for a given burst to be the reported error radius 



divided by 1.58 after Alexandreas et al. |61], but did not employ this optimized scaling 



based on expected numbers of events to then enlarge the already substantial search bin 



radii equal to the quadrature sum of the angular cuts listed in Table |0| and the satellite 
resolution: 



2 satellite position error 



-^search = W ^cnt + ( ^ 



^In the first analysis of tlie data |60|, I cliose A'dirB > 4, '3/ < 10°, wliicfi is sligfitly inferior. 
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Only GRBs within the Northern Hemisphere {9 > 90°) were included in the search. Hard 
cuts were made at the error circle and at ^ = 80° and all events falling within an error bin 
were weighted equally. 

Since particle acceleration is expected to persist only during the internal shocks 
associated with the gamma-ray flash, I confined the search to roughly the duration of each 
GRB. I tallied the number of events, no, falling in the angular search bin beginning 100 
milliseconds before and ending 100 milliseconds after the Tgo time. For those bursts for 
which duration information was not available, I searched over the interval [-2,8] seconds 
around the trigger time. 

I estimated background by averaging events in the same detector angular bin (fixed 
6 and <p) during all other intervals which passed the quality analysis. For most bursts the 
off-source observation was many hundreds or thousands of times the r90 interval and the 
error in the number of expected events, //b, is negligible. 

Results and analysis appear in the next chapter. 
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Chapter 6 



Results and conclusions 



If the glove don't fit, you must acquit! 

Johnnie Cochran 



6.1 Search results 



Figure |6.1| shows the results of the search for neutrino emission during the gamma- 
ray flash of 78 BATSE GRBs at the three cut levels. I compare the chance probability of 
occurrence to that expected with no signal, determined by a Monte Carlo simulation of 
several thousand searches in random angular bins. Results are consistent with no signal. 

Employing the prescription in Montanet et al. ||6^ (page 1281) for a Poisson process 
with background, I place a limit on signal events for each burst and for the combined sample. 
The calculation determines the number of signal events, N, such that any random repetition 
of the experiment with the same expected background, hb, would observe more than the 
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Figure 6.1: Chance probabilities, Pci of on-source searches are consistent at all cut levels 
with that expected with no signal (solid histogram). 



observed number of events no and would have ne < no, all with probability 



g-(MB+Af) ^ 



n=0 



n=0 



0.9. 



(6.1) 



Limits are naturally strongest for short bursts near the horizon at the tightest cut level 
where the Q-factor is relatively high (Figures 6.2 through 6.7). Taking the response in 
AMANDA-BIO given by Monte Carlo simulation (Section [4.2.2 ) and the measured live 
time of 68%, I convert the event limits into a neutrino fluence. With the hard cut I find a 
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combined 90% confidence limit 

El^< 3.8 X 10-Vin{l, TeV cm'^, (6.2) 

per average burst during the gamma-ray flash. This limit is orders of magnitude more 
stringent than the limits of 10^^ and ~ 10~^ upward-going muons cm~^ per burst set 
by the Irvine-Michigan-Brookhaven (1MB) |2^] and the Monopole Astrophysics and Cosmic 
Ray Observatory (MACRO) [^] neutrino detectors. 

6.2 Outlook 

The fluence limit in Equation |6.2| is comparable to the prediction for T = 100, 
and careful modeling of the fireball should permit even the current AMANDA to place a 
useful lower bound on the value of T. Predictions of the GRB neutrino flux will improve as 
more is learned about conditions in GRB fireballs with the launch of a new generation of 
sophisticated satellites. Future GRB missions or further analysis of old data can be expected 
to identify more GRBs, in particular the harder variety accompanied by the sub- TeV to 
TeV emission suggestive of particle acceleration. 

At the preferred cut level used in this search, limits improve roughly as the number 
of bursts in the sample and the addition of the 1998 and 1999 data sets will further constrain 
neutrino emission by a factor of 3 to 4. Enlarging the data set also improves the odds of 
observing the peculiar burst with favorable characteristics. The likelihood of GRB neutrino 
detection is certain to increase as AMANDA and detectors like it continue to expand. 
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Figure 6.2: Limits on neutrino fluence during the gamma-ray flash of 78 GRBs in 1997 

versus burst duration at the tightest cut level. The solid line is the combined fluence limit, 
while dashed lines are model predictions. 
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Figure 6.3: Limits on neutrino fluence during the gamma-ray flash of 78 GRBs in 1997 

versus cos(^) at the tightest cut leveL The sohd hne is the combined fluence hmit, while 
dashed hnes are model predictions. 
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Figure 6.4: Fluence limit versus duration at medium cut level. 
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Figure 6.5: Fluence limit versus cos(^) at medium cut level. 
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Figure 6.6: Fluence limit versus duration at the loosest cut level. 
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Figure 6.7: Fluence limit versus cos(^) at the loosest cut level. 
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